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Resumo: Este trabalho apresenta uma revisao bibliografica em Otimizagao Combinatoria e Programagao Matematica sobre os mode- 
los de localizagao de hubs de maior utilidade para os projetos de instalagao de terminais intermodais de transporte no Brasil. Esta pes- 
quisa foi motivada pela recente alteragao na legislagao sobre concessao, arrendamento e autorizagao de instalagoes portuarias, bem 
como pelo cenario de incentivo govemamental para utilizagao do transporte intermodal no Pais. Aborda-se a classificagao dos tipos de 
metodos, os modelos de formulagao matematica dos problemas e as tecnicas de solugao dos mesmos. 

Abstract: This paper presents a survey on Combinatorial Optimization and Mathematical Programming dealing with the hub location 
models most useful to projects of intermodal terminal installation in Brazil. This research was motivated by the latest changes at the 
legislation on port concession, leasing and installation authorization, as well as by the scenery of government incentive to the use of in- 
termodal transport inside the Country. It approaches the classification of method types, the mathematical formulation models of the 
problems and their solution techniques. 


1. INTRODUQAO 

O Decreto da Presidencia da Republica N°. 6620 de 
29 de outubro de 2008 disciplinou a concessao de por- 
tos, bem como o arrendamento e a autorizagao de ins- 
talagoes portuarias maritimas. Tal ato cria isonomia e 
oportunidades de investimentos para a iniciativa pri- 
vada. Neste sentido, uma das primeiras questoes a se- 
rem resolvidas e a localizagao otima de novos empre- 
endimentos. Assim, cabe ao govemo a determinagao 
das areas costeiras em que se pretende permitir tal ins- 
talagao (definigao que obrigatoriamente passa por atu- 
alizagao bienal). E, por parte de cada investidor, deve- 
ra ser determinado o lugar exato onde localizar seu 
terminal portuario. 

Do lado da administragao publica, o problema foi 
resolvido com a elaboragao do Plano Geral de Outor- 
gas para o Subsetor Portuario (ANTAQ, 2010) e sua 
subsequente aprovagao pela Secretaria Especial de 
Portos. Do lado empresarial, cada investidor tera que 
demonstrar a viabilidade tecnica, operacional e eco- 
nomica do seu projeto especifico, bem como seu im- 
pacto concorrencial, o que passa por um estudo para 
localizagao otimizada do terminal portuario. 

No setor ferroviario, embora nao haja um marco le- 
gislative tao novo como o ja referido, os problemas de 
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localizagao tambem nao sao raros, uma vez que o vo- 
lume crescente de cargas transportadas no setor (con- 
forme se verifica em ANTF, 2010) acaba levando a 
novos investimentos em terminais intermodais por 
parte das concessionarias. 

O incentivo ao transporte intermodal gerado pela 
meta do Govemo Federal de reequilibrio da matriz 
modal de movimentagao de cargas no Pais tambem 
contribui para investimentos em infra-estrutura de 
ter mi nais intermodais. 

No modo aquaviario, a recente abertura gera inclu- 
sive uma expectativa de aumento deste tipo de inves- 
timento, haja vista a demanda reprimida anteriormente 
existente. 

Neste contexto, julga-se apropriado um estudo so- 
bre os metodos de localizagao de ter mi nais de trans- 
bordo, a fim de constituir uma base teorica para anali- 
ses de viabilidade de projetos como os ja referidos. 

Embora a motivagao inicial tenha levado o estudo a 
enfocar sistemas de transportes que por natureza sao 
intermodais, a teoria abordada a seguir tambem pode 
ser aplicada aos casos de terminais unimodais e ate 
mesmo a outros tipos de redes (e. g., telecomunica- 
goes). 

Dados a estrutura fisica dos sistemas de transportes 
e o objeto de estudo exposto, pode-se afirmar que as 
metodologias de localizagao mais adequadas sao a- 
quelas baseadas na Otimizagao Combinatoria - cujo 
processo de modelagem dispoe de uma ferramenta 
fundamental: o grafo. 

O modelo de localizagao especifico a ser utilizado 
em cada aplicagao depende do tipo de fluxo em rede 
analisado. Por exemplo, o escoamento de commodities 
brasileiras ate os portos (para exportagao) nao apre- 
senta uma demanda bem definida entre os produtores 
e os terminais portuarios, de modo que as cargas de 
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uma origem nao necessariamente tem que seguir sem- 
pre para um mesmo porto, isto e, a cada safra, um 
produtor de soja pode variar o porto de destino de sua 
carga sem que isto modifique significativamente sua 
venda. Portanto, neste caso, nao ha um fluxo previa- 
mente determinado entre cada par origem-destino 
(OD), existem apenas as quantidades que se pretende 
exportar a partir de cada origem. Pois a demanda no 
escoamento de produpao ate os portos e bastante sen- 
sivel as mudanpas na oferta de infra-estrutura de 
transporte. 

No sentido inverso (fluxos das importapoes dos por- 
tos ate seus destinos), o quadro se repete, existem so- 
mente as quantidades a serem recebidas em cada pon- 
to no Pais, de modo que os fluxos podem iniciar em 
quaisquer portos disponiveis. Ate mesmo porque, no 
transporte maritimo de longo curso, nao e comum ha- 
ver diferenpa no valor de ffete quando se variam os 
portos a serem utilizados na costa brasileira, isto e, um 
navio graneleiro trazendo fertilizantes do Egito para o 
Brasil cobraria a mesma tarifa quer o porto de desem- 
barque se localizasse em Recife ou Paranagua. 

Para estes casos em que nao existe uma matriz de 
demanda entre pares OD previamente estabelecida 
podem ser empregados tanto os modelos de localiza- 
pao de facilidades classicos quanto os de fluxos em 
rede - ja bastante abordados na literatura sobre Otimi- 
zapao Combinatoria. 

Porem ha os casos onde uma matriz OD de fluxos 
pre-determinados e um dado de entrada do problema. 
Por exemplo, uma cadeia de suprimentos em que um 
produtor de apo atende demandas anuais regulares e 
bem definidas em diversos centros distribuidores do 
Pais. Se este produtor decidir enviar cargas a todos es- 
tes centros exceto para um e em seguida tentar entre- 
gar esta unica carga nao enviada em algum dos demais 
centros “ja abastecidos”, ele nao tera sucesso. E, se 
uma empresa ferroviaria atende este transporte de apo 
e ainda outras cadeias similares, seu planejamento tera 
que considerar sempre uma matriz de fluxos com ori- 
gem e destino ja pre-estabelecidos. 

Para determinado da posipao de novos ter mi nals 
concentradores em sistemas com demanda bem defi- 
nida entre cada par OD (como o deste ultimo exem- 
plo), e necessario aplicar outra classe (mais complexa) 
de modelos, a dos problemas de local izacjao de hubs. 

Tal classe e relativamente ainda pouco explorada e 
tem amp la aplicapao nos projetos de transporte referi- 
dos inicialmente. Desta constatapao surgiu o objetivo 
do presente trabalho : realizar uma analise do estado da 
arte destes modelos de local izado de hubs dircciona- 
da para os referidos projetos. Sendo incluidos neste 
estudo os seguintes topicos: classificapao, modelos de 
formulapoes matematicas dos problemas e as tecnicas 


de solupao dos mesmos. 

Para ilustrar a utilidade destes modelos citam-se al- 
guns exemplos de seu emprego pratico: 

■ local izado de ter mi n a ls intermodais ao longo 
de uma ferrovia, hidrovia ou linha de cabota- 
gem; 

■ local izado de aeroportos publicos; 

■ local iza^ao de terminals unimodais para conso- 
lidapao e desconsolidapao de LCL ( less than 
container load, carga fracionada). 

2. O MODELO DE REDE HUB-AND-SPOKE 

Uma rede hub-and-spoke e aquela em que os vertices 
concentradores {hubs) se conectam aos nos comuns 
por meio de vias para veiculos de menor capacidade 
(conhecidas como spokes), enquanto a conexao entre 
os proprios terminais concentradores e feita por meios 
de maior capacidade (chamados de eixos troncos ou 
shuttle services, Figura 1), os quais proporcionam as 
maiores economias de escala. Portanto um hub recebe 
a convergencia de diversos fluxos (oriundos dos nos 
comuns) e os direciona para outro hub (via um eixo 
tronco) ou os distribui em seguida dirctamcntc a verti- 
ces comuns da rede (via spokes). No sentido inverso, 
analogamente o hub distribui fluxos oriundos de um 
meio de grande capacidade (e/ou de diversos meios de 
menores capacidades). 

O Problema de Localizagao de Hubs consiste em se 
localizar ter mi nais concentradores em uma rede e 



Figura 1 . Esquema de uma rede hub-and-spoke 

tambem em se atribuir a estes hubs os demais vertices 
(terminais comuns), de modo a proporcionar o rotea- 
mento otimo dos fluxos de transporte entre os diversos 
pares de nos, origens-destinos. Nota-se que nao se de- 
fine aqui qual e o fluxo otimo, ja que ha criterios es- 
pecificos para cada tipo de problema de localizagao de 
hub, tais como minima soma (p-medianas), minimo 
maximo caminho (p-centros e maxima cobertura) etc.. 

3. TIPOS DE PROBLEMAS DE LOCALIZAQAO 
DE HUBS (PLH) 

Como cada tipo de problema e definido por uma com- 
binagao de algumas das caracteristicas basicas dos sis- 
temas hub-and-spoke, entao o primeiro passo e des- 
crever cada uma destas, come9ando por aquelas que 
dizem respeito a topologia da rede do problema: 

■ o tipo de alocapao dos terminais comuns ao hub 
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(sendo unica quando um ter mi nal comum pode 
se conectar a apenas um hub da rede e multipla 
quando inexiste tal restricao); 

■ a permissao ou nao de comunicacao direta entre 
nos comuns e 

■ a conectividade entre hubs (total ou parcial). 

Uma quarta caracteristica e a existencia ou nao de 
restricoes de capacidade em vertices ou arcos, o que 
limita o fluxo maximo no caminho que inclui o ele- 
mento “capacitado”. Assim, o problema e dito capaci- 
tado quando ha restricao em algum elemento e nao- 
capacitado no caso contrario. 

Outra caracteristica e o tipo de objetivo da otimiza- 
Cao: localizacao de p -hub medianas (minima soma), 
localizacao de hubs com custos fixos (minima soma), 
localizacao de p -hub centres (minimo maximo cami- 
nho) e localizacao de hubs para cobertura de rede (mi- 
nimo maximo caminho). Uma observacao interessante 
e que, para quase todos os modelos de localizacao de 
hubs, existem versSes analogas de modelos de locali- 
zacao de facilidades classicos. 

A sexta caracteristica identificada nestes problemas 
e o modo de determinacao do numero de hubs na rede, 
que pode ser exogeno ou endogeno (O’Kelly, 1992). 
No primeiro caso, o numero de hubs e definido arbi- 
trariamente pelo analista por razoes, entao tal numero 
se toma mais uma restricao na formulacao do progra- 
ma de otimizacao, que costuma ser denominado pro- 
blema de localizacao de p -hubs. No segundo caso, o 
numero de hubs e uma variavel do problema sendo de- 
finido pelo proprio processo de otimizacao. 

Tendo em vista que a classificacao dos problemas 
pelo tipo de funcao objetivo tem recebido grande a- 
tencao na literatura cientifica, descrever-se-ao a se- 
guir, sucintamente, as propriedades mais importantes 
de cada um dos quatro tipos basicos principals (con- 
forme Campbell, 1994 e Alumur et al.., 2009). 

O Problema de Localizacao de p -Hub Medianas 
consiste em se determinar os locais numa rede para a 
instalacao de um numero p (pre-determinado) de hubs 
e em se atribuir cada fluxo entre dois vertices (origem- 
destino) a um caminho que passe ou nao por hubs, de 
modo que os fluxos de transportes na rede tenham o 
menor custo total possivel. No inicio das pesquisas 
sobre terminals concentradores, este era o tipo de pro- 
blema mais investigado. Para maiores detalhes, deve- 
se consultar Alumur e Kara (2008), Skorin-Kapov et 
al.. (1996) e O'Kelly etal.. (1996). 

O Problema de Localizacao de Hubs com Custos 
Fixos e muito similar ao de p -Hub Medianas e, em 
termos de formulacao, difere deste basicamente em 
dois aspectos: na computacao do custo total (onde sao 
incluidos os custos fixos dos terminals e nao apenas 
os de transports) e na determinacao do numero de 


hubs que e quase sempre endogena. Atualmente, tem 
sido o problema de localizacao de hubs com maior 
numero de publicacoes. 

O Problema de Localizacao de p -Hub Centres e do 
tipo minimo maximo custo. Campbell (1994) define 
as tres classes mais comuns destes problemas: mini- 
mizacao do maximo custo de transports entre todos os 
pares origem-destino; minimizacao do maximo custo 
de transports entre todos os arcos da rede; e minimi- 
zacao do maximo custo de transports sntrs todos os 
pares origem-hub ou hub- dsstino. O rsfsrido autor 
tambem aprsssntou sxsmplos em que os tres casos se 
aplicam respectivamente: transports de cargas urgen- 
tes ou pereciveis, transports de cargas que dependem 
de apoio no terminal (e. g., refrigeracao), sistemas em 
que e necessario evitar longas viagens para os moto- 
ristas dos veiculos de coleta e distribuicao. No Pro- 
blema de Localizacao de p -Hub Centres, e muito co- 
mum a utilizacao do tempo no lugar do custo. Para 
maiores detalhes, ver tambem Kara e Tansel (2000). 

O Problema de Cobertura com Hubs e analogo aos 
problemas classicos de cobertura, em que um ponto de 
demanda e dito coberto por um ponto de distribuicao 
se aquele esta a menos de uma distancia critica deste 
(ou de um custo critico, ou de um intervalo de tempo 
critico etc.). Na cobertura com hubs, Campbell (1994) 
define tres casos em que um par origem-destino (i, j ) 
esta coberto por um par de hubs ( k , m)\ o custo de 
transporte de i para j por meio dos hubs kern nao ex- 
cede determinado limite; o custo em qualquer arco do 
caminho ijkm nao excede o limite; o custo nos arcos 
origem -hub ou /zwh-destino nao excedem o limite. Es- 
tes tres casos podem ser ilustrados respectivamente 
pelos mesmos tres exemplos apresentados para o pro- 
blema de p-Centros. E importante ressaltar que os 
problemas de cobertura se dividem em dois grupos 
principals: o problema de cobertura da rede, cujo ob- 
jetivo e cobrir com hubs todos os fluxos origem- 
destino, minimizando o custo de instalacao do total de 
hubs necessarios para isto; e o problema de maxima 
cobertura, onde se busca cobrir o maximo de fluxos, 
utilizando um numero pre-estabelecido de hubs. Para 
maiores detalhes deve-se recorrer tambem a Kara e 
Tansel (2003). 

4. EVOLUQAO DAS FORMULAQOES DOS 
PLHCOM CUSTOS FIXOS 

Neste item aborda-se o processo de desenvolvimento 
das formulacoes matematicas para problemas de loca- 
lizacao de hubs com custos fixos (tratados neste traba- 
lho por PL/7-CF e conhecidos na literatura cientifica 
por HLP apenas, ja que, na sigla em Ingles e comum 
suprimir a expressao custos fixos). A atencao especial 
que se dedica ao principal problema de localizacao de 
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hubs (entre os quatro classificados por tipo de fungao 
objetivo - descritos no item anterior) se deve ao fato 
de este ser destacadamente o mais utilizado nas apli- 
ca9des que motivaram esta pesquisa (descritas no ini- 
cio do trabalho), uma vez que os modelos de p -Hub 
Medianas e de p -Hub Centres nao incluem os impor- 
tantes custos fixos de instalacjao dos terminais, en- 
quanto o problema de cobertura com hubs tem rara a- 
plica9ao pratica nos referidos casos de interesse. 

Goldman (1969) e conhecido como o primeiro tra- 
balho publicado sobre o Problema de Localiza9ao de 
Terminais Concentradores. Ele apresenta matemati- 
camente a ideia de otimiza9ao do problema e realiza 
algumas demonstra9oes formais. Entretanto O’Kelly 
(1987) e reconhecido como o introdutor da primeira 
formula9ao matematica canonica do PL H. 

4.1. Dados basicos das formulagdes do PLH 

Para se definir matematicamente tal problema e preci- 
so entender a rede de transportes em estudo como 
sendo representada por um grafo G = (N, A), em que 
cada vertice do conjunto N corresponde a um ponto de 
origem e destino dos fluxos (que pode vir a ser sele- 
cionado para a instala9ao de um hub), enquanto os ar- 
cos do conjunto A constituem os caminhos. Define-se 
H <^N como o subconjunto dos vertices considera- 
dos passiveis de receberem a instala9ao de um novo 
terminal hub. Tal subconjunto e normalmente deter- 
minado arbitrariamente e, conforme sua propria defi- 
ni9ao, pode ser o proprio N. W=\ Wy | e definida co- 
mo a matriz de demanda de fluxos Wy de um no i para 
outre no j qualquer, e tal fluxo pode passar por um ou 
dois nos hub intermediaries, compondo-se um cami- 
nh o i — >k — >/77 — (mais frequentemente representado 
pela sequencia ijkni). Note que, apesar de possivel, em 
geral na Engenharia de Transportes, nao se considera 
a utiliza9ao de mais de dois hubs no trajeto entre um 
par origem-destino, pois, tendo em vista que todos os 
hubs concentram muitos fluxos, raramente e necessa- 
rio modelar explicitamente a passagem por mais de 
dois ter mi nais concentradores. 

O custo de transporte da demanda wy por um cami- 
nho ijkm e dado pela equa9ao 1 : 

C ijkm = W ij ' [l ' C ik + a ' C km + b • C mJ ) ( 1 ) 

em que, 

c ib c km, c mj '■ custos unitarios de transportes referentes a 
cada subtrecho; 

X, a, 5 : fatores de desconto devido a economias de 
escala; 

yk e a variavel binaria que define se um hub sera 
instalado ou nao em k (sendo igual a 1 ou 0) e C k o 
custo desta possivel instala9ao. Nas formula9oes tetra- 
indexadas, x, ikm e definida como a variavel referente a 
ffa9ao de Wy que passa pelos hubs kem. 


4.2. Primeiras formulagoes 

A primeira formula9ao do PL/7-CF foi proposta por 
O’Kelly (1992) para o caso nao-capacitado e com alo- 
ca9ao unica, tendo sido concebida na forma de pro- 
grama9ao inteira quadratica, conforme se verifica a 
seguir. 


Min ^ y x ik ■ c ik ■ (O i +D t )+... 

isN kcN 




Sujeita 


^x ik = 1, Vz e N 

(2b) 

y k -£x lt Z0,VkeN 

(2c) 

x ik e {0,l}, Vz e N,\/k e N 

(2d) 

y k e {0,l}, \/k e N 

(2e) 


(2f) 


(2g) 


As variaveis binarias x ik , definem se um vertice i se- 
ra alocado a um hub k (sendo igual a 1 neste caso e 
igual a zero em caso contrario). y k define se um hub 
sera instalado ou nao em k (sendo respectivamente i- 
gual a 1 ou 0). 

A fiunjao objetivo 2a orienta a busca do minimo 
custo total da rede - somando-se custos de transporte 
e custos fixos (instala9ao de terminais etc.). Cada res- 
tri9ao 2b assegura que cada vertice da rede esteja alo- 
cado a exatamente um no hub, enquanto as restr^des 
do tipo 2c garantem que os nos nao sejam alocados a 
hubs que nao venham a ser instalados. As restri9oes 2f 
e 2g definem os totais produzidos e atraidos por cada 
vertice da rede. 

Campbell (1994) elaborou a primeira formula9ao 
“linear” do PL/7-CF cuja ideia geral ainda vem sendo 
utilizada, apesar de ja contar com aprimoramentos in- 
troduzidos por outros pesquisadores. Aquele autor a- 
presentou formula9oes para os quatro casos basicos do 
PL/7-CF - cada caso caracterizado por uma combina- 
9ao das duas caracteristicas fundamentais: capacidade 
(capacitado ou nao-capacitado) e aloca9ao (unica ou 
multipla), ja detalhadas em 3. 

A primeira destas formula9oes (apresentada a se- 
guir) e a do PL//-CF Nao-Capacitado e com Aloca9ao 
Multipla - conhecido na literatura cientifica como 
UMAHLP ( Uncapacitated Multiple Allocation Hub 
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Location Problem). 

Min H E Z Z C iM ' x iM + E C * • Vk ( 3a ) 

Sujeita a, 

EE x iMm =h\/ieN,VjeN (3b) 

< y k , Vi e TV, V/' e TV, Vk e TV, Vw e TV (3c) 

x ijkm < y m , Vi e TV, V/ £N,Vk £ TV, Vm e TV (3d) 
e {0,l}, Vk £ N (3e) 

0 < x ijkm < 1, Vi e TV, V/ e TV, VL e TV, Vm £ AT (3f) 

As variaveis sao as mesmas apresentas em 4.1 e a 
funpao objetivo 3 a orienta a busca do rninimo custo 
total da rede. As restripoes do tipo 3b asseguram que 
todos os fluxos sejam realizados entre todos os pares 
de nos, enquanto as 3c e 3d garantem que os fluxos 
ocorram apenas por meio de hubs realmente instala- 
dos. E importante notar que a determinapao do nume- 
ro total de term in als concentradores e endogena (ex- 
plicado em 3). 

Para se obter a formulapao de Campbell (1994) da 
versao “Capacitada” do PL//-CF com Aloca9ao Mul- 
ti pi a, basta adicionar a formula9ao 3a-3f a seguinte 
restripao: 

Z2>v • + x ‘m )~ *#*] * r * • y^ k e N (3g) 

em que e a capacidade maxima de um hub localiza- 
do no vertice k. Esta desigualdade garante que, caso 
haja um hub instalado em k, o somatorio de todos os 
fluxos que passam por tal hub nao excedera a sua ca- 
pacidade r t . 

Por outro lado, caso se deseje manter o PL//-CF 
como Nao-Capacitado porem alterar o tipo de aloca- 
pao (transformando-o no tipo Alocapao Unica), sao 
necessarias duas alterapoes na formulapao 3a-3f: a 
substituipao das restripoes 3c e 3d pela restripao 3h e a 
inserpao das restripoes 3i e 3j (apresentadas a seguir). 
A variavel binaria z ik determina se o vertice i esta alo- 
cado ao hub k (sendo igual a 1 neste caso e zero em 
caso contrario). 

z ik < y k , Vi e N, Vk e TV (3h) 



Vi e N,Vk e TV 


z ik e {0,l}, Vi s N,Vk & N (3j) 


A restripao 3h garante que um vertice i somente po- 
de ser alocado a um hub em k caso este ultimo seja 
instalado. E a desigualdade 3i obriga que todos os flu- 
xos que cheguem ou saiam de i o fapam pelo hub em 
k, ao qual i esta alocado. 

Finalizando as formulapoes de Campbell (1994) pa- 
ra o PL/V-CF, descreve-se aqui a modificapao necessa- 
ria para, a partir do terceiro caso exposto (PL/V-CF 
“Nao-Capacitado” e com Alocapao Unica), chegar ao 
PL//-CF “Capacitado” e com Alocapao Unica. Para 
tanto, basta que se adicione ao problema nao- 
capacitado a restripao a seguir, a qual garante que, ca- 
so um vertice i esteja alocado a um hub instalado em 
k, o somatorio de todos os fluxos que chegam e saem 
de i nao excedera a capacidade daquele hub, T k . 

EEk + * T k ■ y k ,Vk £ TV (3k) 

4.3. Modificagoes propostas para as 
primeiras formulagoes 

Tendo sido abordadas as formul apdes originais dos 
quatro casos principais do Problema de Local izapao 
de Hubs com Custos Fixos (PL/V-CF), a partir deste 
ponto, o texto sera dedicado ao estudo da evolupao 
das formulapoes para os casos de alocapao “multipla” 
(definida no item 3), haja vista que os problemas pra- 
ticos que motivaram esta pesquisa estao relacionados 
a local izapao de terminals intermodais que quase sem- 
pre tem o modo rodoviario como a altemativa princi- 
pal (se nao unica) de distribuipao e coleta de cargas 
para os hubs. E os operadores desta modalidade mode- 
lo trabalham com total liberdade e flexibilidade na es- 
colha de suas rotas, de tal modo que a aloca9ao unica 
nao retrataria com fidelidade a realidade destes siste- 
mas. 

Na alocapao multipla, por exemplo, um no i, que 
envia cargas a dois nos distintos j e /, pode utilizar um 
par de hubs km como intermediario para o trajeto ij e 
outro par pq para o trajeto il. Ja no caso de alocapao 
simples, cada no i da rede e cativo de um unico termi- 
nal concentrador k (normalmente o que apresenta o 
menor custo de acesso ao mesmo). Vale tambem res- 
saltar que a alocapao unica nao necessariamente mi- 
nimiza os custos totais de uma rede em um problema 
de local izapao de hubs, mas a alocapao multipla sim. 

Skorin-Kapov et al.. (1996), estudando o Problema 
de p -Hub Medianas, descobriram desigualdades mais 
fortes ( tighter ) que as restripoes dos tipos 3c e 3d de 
Campbell (1994). Apesar de se tratar de um problema 
um pouco diferente do PL//-CF, as melhorias obtidas 
tambem sao aplicaveis a este ultimo. Assim, aqueles 
autores substituiram as referidas restripoes (conforme 
Tabela 1), obtendo uma nova formulapao mais aperta- 
da (restringida), com a vantagem de que esta modifi- 
capao nao implicou em um aumento do numero de va- 
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Tabela 1 . Substitutes de restrigoes propostas por Skorin-Kapov, Hamacher etal. 


Restrigao Campbell (1994) Skorin-Kapov et al. (1996) Hamacher et al. (2004) 

3c X ijbn~yk 'Yj X ijkm~yk 

Z X ijkm + L Xi M — y* 

3d x ijkm s y m E*^ - y m 


riaveis no problema, e ainda reduziu o numero de res- 
trigoes para (2n 4 - 2n 3 ). 

Hamacher et al. . (2004), a partir de uma analise po- 
liedrica do problema classico de Localizagao de Faci- 
lidades (PLF) Nao-Capacitado, descobriram facetas 
do espa 90 de solugao do PL// Nao-Capacitado e, tam- 
bem, as dimensoes destas. O inverso tambem foi con- 
seguido, isto e, determinaram facetas para o PLF a 
partir do estudo do Problema de Localizagao de Hubs. 
Esta inovagao propiciou mais uma modificagao nas 
restrigoes 3c e 3d, permitindo a fusao das mesmas em 
uma unica (conforme Tabela 1) e tomando a formula- 
gao ainda mais apertada (tighter). 

Marin et al.. (2006) criticam a formulagao 3a-3f 
(mesmo alterada por Skorin-Kapov et al.., 1996). A- 
firmam que esta formulagao “nao” garante que uma 
rota (definida por quatro indices, ijkm - origem, hubl, 
hub 2 e destino) possuira no maximo dois hubs (limite 
usual nestes problemas). Eles desenvolveram uma no- 
va formulagao (para o PL//-CF Nao-Capacitado e com 
Multipla Alocagao) composta de restrigoes do tipo fa- 
cetas do poliedro de solugoes do problema. Utiliza- 
ram, para isso, uma tecnica criada por Padberg (1973) 
apud Marin et al.. (2006) baseada em restrigoes de 
cobertura de conjunto (set packing constraints ), nas 
quais a soma de uma serie de variaveis e menor ou i- 
gual a 1 . Contaram tambem com o fato de o problema 
ser do tipo Nao-Capacitado, onde as variaveis x, jkm sao 
binarias. Assim, obtiveram a formulagao final 4a-4j. 


Min H ZEE c iM ■ x iM +Z C * a* 

ieN jeN ksN meN keN 

(4a) 

EEv =1 ’ Vi gN,\/J gN 

keN meN 

(4b) 

E H X ‘M + y> < l ’ e N > vj e N ’t * j 

keNmeN 

(4c) 

E + yj - l ’ yi e N > v -/‘ e N ’i * j 

keNmeN 

(4d) 

h+ E 

(4e) 


+ E*iM * y*> Vi e N . V/ e N, Vk e N (4f) 

x kkbn = x kbnk = o, VA: e N, Vw e N, m * k (4g) 
x kjmk = x jkkm = 0, V* e N, Vw e N, V/ e N, 

k,,nk jkkm ( 4h ) 

m*k A j *k 

y k m {0,l}, VA: e N (4i) 

0 < x ljhn , Vi e N, V/ gN ,Vk&N,Vm&N (4j) 

Esta formulagao apresenta somente (n 3 + 3n 2 - n) 
restrigoes (sendo n o numero de vertices da rede), o 
que e bastante inferior a formulagao original de 
Campbell (1994), mesmo melhorada por Skorin- 

Kapov et al.. (1996), que ainda possui (2n 4 - 2n 3 ) res- 
trigoes. 

Segundo Alumur e Kara (2008), com esta formula- 
gao 4a-4j, Marin et al.. (2006) generalizaram as for- 
mulagdes existentes do PLH-CF Nao-Capacitado e 
com Multipla Alocagao e relaxaram a restrigao de que 
a estrutura de custos da rede tenha que satisfazer a de- 
sigualdade triangular. Alem disso, utilizando algumas 
tecnicas de pre-processamento (descritas no item 

3.5.5), eles produziram uma formulagao mais forte 
(mais apertada, tighter) que superou em desempenho 
todas as anteriormente publicadas. 

4.4. Outros tipos de formulagoes 

As formulagoes apresentadas em 4.2 e 4.3 utilizam va- 
riaveis bi-indexadas ou tetra-indexadas para a defini- 
gao do roteamento dos fluxos (por exemplo, Xy ou 
Xy/an). Existem tambem, embora em menor numero, as 
formulagoes com variaveis de alocagao tri-indexadas 
(Ximj) e outras com variaveis baseadas nos arcos da re- 
de (ao inves dos vertices). Alem disso, ha outras for- 
mulagoes que introduzem importantes contribuigoes 
mas que, em termos de forma matematica, apresentam 
apenas pequenas modificagoes em relagao as formula- 
goes descritas em 4.2 e 4.3. Esta segao aborda estes 
casos. 

Ebery et al. . (2000) e Boland et al. . (2004) elabora- 
ram formulagoes para o PL//-CF com Alocagao Mul- 
tipla (tanto Capacitado quanto Nao-Capacitado), utili- 
zando variaveis tri-indexadas e obtiveram uma consi- 
deravel redugao do numero de variaveis e restrigoes, 
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chegando a respectivamente (2n 3 + n) e (2n 2 + 2n). 
Cabe ressaltar que o arcabou90 matematico da maior 
parte desta inovasao foi concebido em Ernst e Krish- 
namoorthy (1998). 

Marin (2005) e Marin et al. . (2006) afirmam que es- 
ta referida formula^o possui as vantagens do menor 
tamanho e do tempo de sokujao muito mais curto em 
rela9ao as formula9des tetra-indexadas (como as 4a-4j 
e anteriores), porem apresentam a desvantagem de os 
limites encontrados na relaxa9ao linear serem mais 
folgados, isto e, piores. Marin et al.. (2006) destacam 
tambem que tal formula9ao falha quando a configura- 
pao da rede nao satisfaz a desigualdade triangular, ja 
que, neste caso, a formula9ao nao garante que cada ro- 
ta entre origem e destino possuira no maximo dois 
hubs (limite usual nestes problemas) e, alem disso, 
pode permitir que, em algumas rotas, certos vertices 
nao -hubs sejam atravessados. Para resolver tais ques- 
toes, Marin et al. (2006) propuseram a inclusao de ou- 
tras restr^des, conforme tambem o fizeram para as 
formula9oes tetra-indexadas do item 4.3. 

Marin (2005) desenvolveu uma formula9ao tri- 
indexada especifica para um PL/7-CF-AM Capacitado 
que permite o particionamento de cada fluxo entre 
dois vertices (wy) em diversas rotas. Este problema foi 
denominado pelo autor como Splittable Capacitated 
Multiple Allocation Hub Location Problem. 

Rodriguez-Martin e Salazar-Gonzalez (2008) estu- 
daram um problema de local izapao de hubs com cus- 
tos fixos e capacitado; ao qual deram o nome Capaci- 
tated Hub Location Problem ( CHLP ), cuja formula- 
9ao tambem e bastante diferente das formula9oes tetra, 
tri e bi-indexadas referidas, ja que estas ultimas se ba- 
seiam nas quantidades de fluxos em cada “rota” (re- 
presentadas por variaveis do tipo x ijkm ) enquanto a dos 
referidos autores se baseia apenas nas quantidades de 
fluxos em cada “arco” da rede (representadas por va- 
riaveis do tipo x a ). 

Yang (2009) apresenta uma importante contribui9ao 
para modelar com maior precisao os PL/7-CF-MA nos 
casos em que a demanda de transporte e caracterizada 
por consideravel variapao sazonal, como no transporte 
aereo. Assim, embora a formula9ao esteja baseada nas 
ja referidas variaveis tetra-indexadas, ela se diferencia 
por constituir uma programapao estocastica e multi- 
estagio. 

Alumur et al. . (2009) propuseram a mudan9a de um 
paradigma dos modelos de localiza9ao de hubs em re- 
des, onde sempre se considerava que, tendo sido insta- 
lado dois novos hubs em uma rede, uma via de liga9ao 
entre os mesmos tambem obrigatoriamente estaria ins- 
talada (ver Figura 1). Portanto Alumur et al.. (2009) 
incluiram, em cada fmujao objetivo de problema estu- 
dado, o custo de instala9ao desta via. 

Camargo et al.. (2009) modificaram a formula9ao 


para o PL//-CF-MA Nao-Capacitado proposta por 
Hamacher et al.. (2004, referida no item 4.3). Tal mo- 
difica9ao tomou possivel considerar no referido pro- 
blema os custos de congestionamento de fluxos em 
cada hub. 

Correia et al.. (2010) elaboraram uma extensao do 
Problema de Localiza9ao de Hubs Capacitado e com 
Aloca9ao Unica em que a principal inova9ao e o fato 
de o limite de capacidade de cada hub ser tambem 
uma incognita do problema. 

Ishfaq e Sox (2010) apresentaram um modelo do 
problema de p -hub medianas “com custos fixos” bas- 
tante rico em detalhes para o emprego no transporte 
intermodal. Como a unica diferenpa deste problema 
para o PL//-CF e a restr^ao que determina o numero 
exato de hubs a ser instalado, vale a pena tratar deste 
caso tambem pela sua grande contribui9ao para as re- 
des intermodais. A formulapao permite uma modela- 
gem mais precisa, robusta e melhor que as formula- 
9oes existentes para o caso da aplica9ao em redes lo- 
gisticas intermodais. Apresenta parametros e variaveis 
bastante detalhados de modo a considerar diferentes 
tipos de rotas (diretas, via hubs com aloca9ao multipla 
e via hubs com aloca9ao unica), diferentes modos de 
transporte (rodoviario, ferroviario e aereo), custos fi- 
xos (incomuns nos p -hubs medianas), custos de trans- 
bordo intermodal, economia-de-escala nao linear e 
restri9oes de tempo de servi90. Com tantos detalhes, 
alguns tipos de variaveis chegam a possuir cinco indi- 
ces. 

5. MODELOS DE SOLUQAO 

Nesta se9ao apresentam-se, sucintamente, os metodos 
mais eficientes encontrados na literatura cientifica e 
destinados a solu9ao do Problema de Localiza9ao de 
Hubs com Custos Fixos. Procurou-se esgotar o assun- 
to para o caso de alocapao multipla, mas nao para o de 
aloca9ao unica (conforme justificado em 4.3), embora 
estejam descritos tambem modelos de sokujao do 
PL//-CF de Aloca9ao Unica. 

No inicio das pesquisas, era nitida a preferencia por 
heuristicas, o que se deve provavelmente a alta com- 
plexidade do problema (NP-Completo) e a existencia 
de pouco conhecimento sobre o mesmo na epoca. A- 
tualmente, conforme mostra Camargo et al. . (2008), a 
pesquisa por metodos exatos e mais ffequente e tem 
base no argumento destes autores de que o planeja- 
mento e projeto de redes hub-and-spoke cria decisoes 
sobre a aplica9ao de grandes quantidades de recursos, 
as quais tem consideravel impacto nos custos totais 
destas redes. O que justifica o gasto de algumas horas 
quando se trata de obter a sokujao otima para um pro- 
blema de localiza9ao de hubs. 

O’Kelly (1992), o primeiro artigo sobre o problema 
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“com custos fixos”, apresenta um procedimento com- 
posto por tres passos principals: uma rotina para de- 
terminapao de um limite superior (uma solupao viavel 
nao necessariamente otima), um segundo procedimen- 
to para calculo de um limite inferior justo e um passo 
final de melhoria do limite superior com base no infe- 
rior. 

Aykin (1994), para o PLiZ-CF-AM Capacitado, cri- 
ou uma heuristica e um algoritmo de branch-and- 
bound e, a partir destes, obteve os limites inferiores 
por relaxapao lagrangeana com otimizapao de subgra- 
dientes. Aykin (1995), para o PL/7-CF-AM Nao- 
Capacitado, tambem utilizou um algoritmo de branch- 
and-bound e o combinou com uma heuristica gulosa 
baseada na meta-heuristica Simulated Annealing. 

Klincewicz (1996) propos um modelo baseado no 
Metodo Dual Ascendente que pode ser entendido co- 
mo uma analogia (para hubs) do algoritmo de Erlen- 
kotter (1978). Klincewicz compos seu metodo em tres 
partes: um algoritmo dual ascendente, um ajuste dual 
e um algoritmo de branch-and-bound. O artigo apre- 
senta resultados de testes computacionais com redes 
de ate 25 vertices, que chegaram a solupdes otimas ou 
proximas delas (95%). 

Mayer e Wagner (2002) desenvolveram importantes 
modificapoes para a metodologia de Klincewicz 
(1996) e apresentaram um artigo mais elucidativo que 
o original. O algoritmo modificado tambem apresen- 
tou desempenho bastante superior ao de Klincewicz. 
Entretanto, ao compararem o HubLocator com o paco- 
te comercial CPLEX, nao foi possivel confirmar a su- 
perioridade de algum dos dois. 

Sung e Jin (2001) tambem utilizaram o Metodo Du- 
al Ascendente para resolver um problema de localiza- 
pao de hub nao-capacitado com a particularidade de 
que a rede deveria ser dividida em clusters (conjuntos 
fechados de vertices) e que houvesse exatamente um 
hub em cada cluster. Sendo assim, o problema, por 
natureza, ja toma a forma de Alocapao Unica. Outro 
detalhe importante e que eles consideraram a possibi- 
lidade de ocorrer fluxos diretos entre origem e destino 
sem passar por hubs. Este metodo de Sung e Jin 
(2001) e composto por um procedimento dual ascen- 
dente e outro de ajustamento dual. 

Canovas et al. . (2007) tambem criaram uma heuris- 
tica baseada no Metodo Dual Ascendente e propuse- 
ram um algoritmo original para resolver o PL/7-CF- 
AM Nao-Capacitado. Antes das iterapoes do algorit- 
mo base, e realizado um pre-processamento que, se- 
gundo os autores, melhora bastante a eficiencia do 
metodo. O algoritmo que coordena todas as funpoes e 
do tipo branch-and-bound e, em cada no desta enume- 
rapao implicita, uma rotina heuristica dual ascendente 
e executada como a ferramenta fundamental do pro- 
cesso. Uma observapao interessante dos autores e que 


as instancias de problemas mais dificeis sao aquelas 
em que ha simetria nos custos de transporte e baixos 
valores de a (coeficiente de desconto para custos inte- 
rhub). 

E importante destacar que, ate o momento, entre to- 
dos os metodos publicados, este metodo de Canovas 
et al. . (2007) e o segundo com o melhor desempenho 
computacional para o PL7/-CF-AM, resolvendo ins- 
tancias de ate 120 vertices em cerca de 4.700 segun- 
dos para problemas assimetricos e 20.600 segundos 
para casos simetricos. O metodo que o supera em efi- 
cacia e o de Camargo et al.. (2008), que sera descrito 
adiante. 

Aqui, o texto volta a ordem cronologica das publi- 
capoes, que foi interrompida porque os tres ultimos ar- 
tigos comentados tem relai^ao direta de tema com o de 
Klincewicz (1996, Metodo Dual Ascendente). 

Ernst e Krishnamoorthy (1999) resolveram um mo- 
delo do PL7/-CF Capacitado e com Alocapao Unica. 
Eles propuseram um metodo gerenciado por um algo- 
ritmo de branch-and-bound que utiliza a relaxapao li- 
near como base das iterapoes e uma combinapao das 
meta-heuristicas Simulated Annealing e Random Des- 
cent para obter limites superiores e buscar a solupao 
de forma mais rapida. Os autores tambem criaram tres 
tecnicas especificas de pre-processamento da formula- 
pao. 

Ebery et al.. (2000), para a solupao do PL/7-CF-AM 
Capacitado, criaram uma heuristica que utiliza um efi- 
ciente algoritmo (baseado na tecnica do caminho mais 
curto - Shortest Path) de Ernst e Krishnamoorthy 
(1998) para gerar solupdes iniciais. Como tal algorit- 
mo foi desenvolvido para problemas nao-capacitados, 
ele acaba gerando solupoes que nao atendem as restri- 
poes de capacidade do problema de Ebery et al.. 
(2000). Assim, estes ultimos autores incluiram uma 
rotina para, a cada solupao inicial, redistribuir os flu- 
xos que estiverem violando restripoes de capacidade. 
Entao, havendo uma solupao inicial valida, e aplicada 
uma rotina gulosa de busca na vizinhanpa desta solu- 
pao para encontrar melhoras, isto e, uma rotina tenta 
iterativamente trocar a local izagao de algum hub da 
solupao corrente, verificando a melhora da funpao ob- 
jetivo. Em caso de melhora, a nova solupao encontra- 
da passa a ser a inicial e uma nova iterapao e proces- 
sada, ate que se atinja algum criterio de parada. 

Wagner, B. (2007) estudou o mesmo problema de 
Sung e Jin (2001) - o Problema de Local izagao de 
Hubs Nao-Capacitado com Clusters (consequente- 
mente, de Alocapao Unica). O metodo exato que pro- 
pos superou em desempenho a heuristica dual ascen- 
dente de Sung e Jin, nao somente no aspecto de pro- 
duzir solupoes exatas mas principalmente por utilizar 
um tempo de processamento significativamente me- 
nor. 
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Camargo et al. (2008) criaram o algoritmo que, atu- 
almente, apresenta o melhor desempenho na solupao 
“exata” do PL/f-CF-AM Nao-Capacitado. O procedi- 
mento e baseado na classica Decomposipao de Ben- 
ders e e capaz de resolver instancias de ate 200 verti- 
ces em menos de 10.000 segundos de processamento. 
Os autores tambem avaliaram mais duas versoes do 
algoritmo de Benders (nao elaboradas pelo proprio 
Benders). Os resultados indicaram que o procedimen- 
to classico superou os demais, especialmente nas ins- 
tancias de maior porte. 

Outro algoritmo baseado no Metodo de Decompo- 
sipao de Benders foi proposto em Camargo et al. 
(2009) para a solupao do PL/7-CF-AM Nao- 
Capacitado sob Congestionamento (problema descrito 
em 4.4). Embora o problema seja uma complexa vari- 
ante nao-linear do PL//-CF-AM, o algoritmo resolveu 
instancias de ate 81 vertices e mostrou desempenho 
superior ao de um software lider comercial em pro- 
grama?ao inteira nao-linear. 

Rodriguez-Martin e Salazar-Gonzalez (2008) testa- 
ram quatro metodos distintos para resolver o Capaci- 
tated Hub Location Problem ( CHLP ), tendo encontra- 
do melhor desempenho em uma heuristica baseada na 
combinapao da meta-heuristica VNS ( Variable Neigh- 
borhood Search) com a relaxapao linear do problema 
(que foi capaz de tratar instancias grandes e fora do 
alcance dos tres metodos anteriores, apesar de nem 
sempre atingir a solupao otima). 

Tres trabalhos bastante recentes publicaram meto- 
dos de solupao para o PL//-CF com Alocacjao “Uni- 
ca”. Silva e Cunha (2009) apresentam modelos de so- 
lu9ao para o caso nao-capacitado deste problema. Tres 
dos modelos se baseiam na heuristica Busca-Tabu 
Multistart, enquanto o quarto e baseado na Busca- 
Tabu Integrada de dois estagios. Tratam-se de meto- 
dos simples com facilidade de implementapao e ao 
mesmo tempo robustos na qualidade das solupoes en- 
contradas segundo afirmam os autores quando repor- 
tam seus experimentos com alguns conjuntos de dados 
intemacionalmente conhecidos. 

O segundo dos tres trabalhos aborda o caso “capaci- 
tado” do PL//-CF com Alocapao Unica. Randall 
(2008) desenvolveu quatro modelos de solupao basea- 
do numa heuristica do tipo Colonia de Formigas com 
eficientes estrategias de busca local e pre- 
processamento. O referido autor afirmou ter encontra- 
do resultados de qualidade tambem com um dos con- 
juntos de dados utilizados por Silva e Cunha (2009). 

Contreras et al. (2009) - no terceiro dos referidos 
tres trabalhos - elaboraram tambem para o caso capa- 
citado um metodo em que a tecnica de Relaxapao La- 
grangeana e especificamente direcionada para a estru- 
tura do problema em questao, decompondo-o em pro- 
blemas menores e gerando um intervalo de dualidade 


(entre os limites superior e inferior da solupao) que 
posteriormente e bastante diminuido por uma heuristi- 
ca especifica baseada em testes de redupao. Tal que os 
intervalos finais de dualidade de todas as instancias 
testadas pelos autores nao superaram 3,4%. 

6. MODELAGEM DA ECONOMIA DE ESCALA 

Desde que nao haja congestionamentos significativos 
na rede, quanto maior for a quantidade de fluxos num 
determinado arco da rede, menor tende a ser o custo 
unitario de transporte neste arco. Este conceito de e- 
conomia vem sendo aplicado na modelagem de redes 
com ter mi nals concentradores, por meio da inclusao 
de fatores de desconto nos custos de transporte de ca- 
da arco (conforme a equapao 1). 

A maioria dos trabalhos sobre PL H existentes con- 
sidera estes fatores de desconto como constantes do 
tipo 0<a<l. Contudo algumas pesquisas ja melhora- 
ram tal modelagem, aproximando mais este fator da 
realidade, pois o consideram uma funpao da quantida- 
de de fluxo. 

O’Kelly e Bryan (1998), modelando o PL/7-CF- 
AM, propuseram, para o calculo do custo de transpor- 
te em um arco km, uma funpao nao-linear concava 
(semelhante a CT NL km da Figura 2), que aumenta o va- 
lor do custo a uma taxa decrescente a medida que os 
fluxos crescem. As demais funpoes da Figura (CT SE km 
e Cf' km ) representam respectivamente o referido custo 
total sem economia de escala e com economia aplica- 
da por um fator a constante (referido anteriormente). 
Flomer e O’Kelly (2001) tambem propuseram outra 
funpao de custo nao-linear. 

Racunica e Wynter (2005) apresentaram outro mo- 
de lo para o PL/f-CF-AM Nao-Capacitado no trans- 
porte intermodal de cargas, onde utilizaram uma fun- 
<?ao de custo nao-linear concava nos trechos interhubs 
e nas ligapoes entre hub e destino (excetuando apenas 
as ligapoes entre origem e hub). Essa fungao e a 
CT W km da Figura 2. Assim, para resolver o problema 
de forma linear, eles fizeram uma aproximagao de ca- 
da funpao concava por meio de diversas funpoes linea- 
res e identificaram algumas propriedades poliedricas 
do modelo linearizado. Entao desenvolveram duas 
heuristicas de redupao de variaveis para a solupao do 
problema e apresentaram um estudo de caso sobre 
uma rede europeia transalpina de transporte de cargas. 

7. CONSIDERAQOES FINAIS 

Foi realizada, na literature cientifica sobre Otimizagao 
Combinatoria e Programapao Matematica, uma revi- 
sao bibliografica sobre os modelos de local izagao de 
ter mi nals concentradores (hubs). Esta pesquisa foi di- 
recionada para o objetivo de aplicar a referida teoria 
em projetos de instalapao de terminais intermodais de 
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transporte do Brasil, contribuindo pro-ativamente para 
a demanda de planejamentos e estudos de viabilidade 
a serem provavelmente requeridos apos a recente pu- 
blicaijao do Decreto 6620 da Presidencia da Republi- 
ca, que passa a vigorar ja num cenario de incentivo 
govemamental para utilizaqao do transporte intermo- 
dal no Pais. 

A pesquisa abordou principalmente a classificaqao 
dos tipos de metodos, seus modelos de formulapao 
matematica e as tecnicas de soluijao dos mesmos; ten- 
do mostrado a diversidade de solupdes para problemas 
desta natureza, bem como a existencia tambem de bas- 
tante oportunidade de adicional desenvolvimento nes- 
ta area (pouco explorada em relapao a teoria classica 
de Localizapao de Facilidades). O trabalho tambem 
leva a conclusao de que os metodos estudados sao 
muito adequados ao planejamento dos sistemas de 
transportes objeto do estudo. 
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